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ABSTRACT 

Jet driving and fragmentation process in collapsing primordial cloud are stud- 
ied using three-dimensional MHD nested grid simulations. Starting from a ro- 
tating magnetized spherical cloud with the number density of n c ~ 10 3 cm -3 , we 
follow the evolution of the cloud up to the stellar density n c ~ 10 22 cm -3 . We cal- 
culate 36 models parameterizing the initial magnetic and rotational energies (70, 
Po). In the collapsing primordial clouds, the cloud evolutions are characterized 
by the ratio of the initial rotational to magnetic energy, 7o//?o- The Lorentz force 
significantly affects the cloud evolution when 70 > /3o, while the centrifugal force 
is more dominant than the Lorentz force when /3 > 7o- When the cloud rotates 
rapidly with angular velocity of fi > 10~ 17 (?2 c /i0 3 cm" 3 ) 2 / 3 s -1 and /3 > 70, 
fragmentation occurs before the protostar is formed, but no jet appears after the 
protostar formation. On the other hand, a strong jet appears after the proto- 
star formation without fragmentation when the initial cloud has the magnetic 
field of B > 10" 9 (n c /10 3 cm" 3 ) 2 / 3 G and 70 > (3 . Our results indicate that 
proto-Population III stars frequently show fragmentation and protostellar jet. 
Population III stars are therefore born as binary or multiple stellar systems, and 
they can drive strong jets, which disturb the interstellar medium significantly, as 
well as in the present-day star formation, and thus they may induce the formation 
of next generation stars. 

Subject headings: binaries: general — cosmology: theory — early universe — ISM: 
jets and outflows — MHD — stars: formation 



1. INTRODUCTION 

Magnetic field plays an important role in present-day star formation. Observations in- 
dicate that molecular clouds have the magnetic field strengths of order ~ fiG and magnetic 
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energies comparable to the gravitational energies (jCrutcherlll999l ). These strong fields signif- 
icantly affect star formation processes. For example, protostellar jets, which are ubiquitous 
phenomena in the star-forming region, is considered to be driven f rom the protostar by the 
Lorentz force (IBlandford fc Payne! Il982l ; iPudritz &: Norman! Il986l ). The jet influences the 
gas accretion onto the protostar and disturbs the ambient medium. In addition, t he angular 
mome ntum of the cloud is removed by the magnetic braking and protostellar jet. iTomisaka 
( 120001 ) showed that, in his two dimensional MHD calculation, 99.9% of the angular momen- 
tum is transferred from the center of cloud by the magnetic field. This removal process 
of angular momentum makes the protostar formation possible in the parent cloud that has 
much larger specific angular momentum than that of the protostar. On the other hand, so 
far, magnetic effects in primordial gas clouds have been ignored in many studies because 
the magnetic field in the early universe is supposed to be extremely weak. However, recent 
studies indi c ate th at moderate strength of magnetic fields exists even in the early universe. 
Ichiki et al.l (120061 ) showed that cosmological fluctuations produce magnetic fields before the 
epoch of recombinat ion. These fields are large enough to seed the magnetic fields in galaxies. 
Langer et al.l (120031 ) proposed the generation mechanism of magnetic fields at the epoch of 
the reionization of the universe. They found that magnetic fields in the intergalactic matter 
are amplified up to ~ 10~ n G. These fields can therefore increase up to ~ 10~ 7 — 10~ 8 G 
in the first collapsed object having the number density n ~ 10 3 cm~ 3 . These fields may 
influence the evolution of primordial gas clouds and the formation of Population III stars. 

Under the spherical symmetry including the hydrodynamical radiative transfer, the 
star form ation process are carefully investiga t ed by many authors both in the present- 
day (e.g.. iLarsonl Il969l ; iMasunaga fc Inutsukal 120001 ) and primordial star formation (e.g., 



Omukai fc Nishi 19981 ). A significant difference between present-day and primordial star for- 
mation exists in the thermal evolution of the collapsing gas cloud because of difference in 
the abundance of dust grains and metals. In present-day star formation, gas temperature in 
molecular clouds is ~ 10 K. These clouds collapse isothermally with polytropic index 7 ~ 1 
(P oc p 7 ) in n c < 10 11 cm -3 , then the gas becomes adiabatic (7 ~ 7/5) at n c ~ 10 11 cm -3 , 
and an adiabatic core (or the first core) is formed, where n c denotes the central number 
density of the collapsing cloud. When the number density reaches n c ~ 10 16 cm -3 , the 
molecular hydrogen is dissociated (7 ~ 1.1), and the cloud begins to collapse rapidly again. 
When the number density reaches n c ~ 10 20 cm -3 , the equation of state becomes hard again 
(7 ~ 5/3). The protostar is formed at n c ~ 10 21 cm~ 3 . On th e other hand, primordial gas 
clouds have temperature of ~ 200 — 300 K at n c — 10 3 cm" 3 (lOmukail l2000l ; lOmukai et al. 



20051 : iBromm et al.ll2002l : lAbel et al.ll2002l ; lYoshida et al.l 120061 ) . These clouds collapse keep- 
ing 7 ~ 1.1 for a long range of 10 4 < n c < 10 16 cm -3 . After the central density reaches 



Tic 



10 16 cm 3 , the thermal evolution in the primordial collapsing cloud begins to coincide 
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with that in present-day cloud (for detail, see Fig. 1 of lOmukai et al.ll2005l ). The differ- 
ence in thermal evolution between present-day and primordial clouds arises in the period of 



11 r 



< 10 cm" . In present-day star formation, two distinct flows (molecular outflow and op- 
tical jet), which are frequently observed in star forming regions, are driven by the respective 
core; the o utflow is driven by the adiabatic cores (the first core) and the jet is driv en by the 
protostar jTomisakalbooilMachida et alil2005bl . l2007bl : Eaneriee fc Pudritzll2006h . In addi- 
tion, many numerical simulations have shown that fragmentation of the adiabatic core, and 
thus possible formation of binary or multiple stellar systems. In contrast, primordial proto- 
stars are formed directly without the prior formation of adiabatic core, because the thermal 
pressure smoothly increases keeping 7 ~ 1.1 for n c < 10 16 cm -3 . Even when the primordial 
cloud is strongly magnetized, the wide-angle outflow which corresponds to the present-day 
molecular outflow (that is driven from the adiabatic core formed at n c ~ 10 11 cm -3 ) may 
not appear, but the well-collimated jet (that is driven from the protostar) is ejected at 
n c > 10 21 . In addition, it is expected that, in the collapsing primordial cloud, fragmen- 
tation rarely occurs for n c < 10 16 cm~ 3 , because the cloud continues to collapse and the 
perturbation inducing fragmentation cannot grow sufficiently 

Another major difference between present-day and primordial star formation exists in 
the magnetic evolution. In present-day star formation, the neutral gas is coupled well with 
the ions in n c < 10 12 cm -3 and n c > 10 15 cm -3 , while the magnetic fi eld are dissipated b y 
the Ohmic diss i pation in 10 12 cm -3 ; 
Machida et al.l (l2007al ) showed that 



n c < 10 5 cm 3 (for details, see iNakano et al.l 120021 ) . 



99% of the magnetic field in the collapsing cloud is 
dissipated in 10 12 cm -3 < n c < 10 15 cm~ 3 . On t he other hand , the magne tic evolution in a 
primordial gas cloud is carefully investigated by lMaki &: Susal (12004 120071 ). and they found 
that the magnetic field strongly couples with the primordial gas during all the c ollapse phase, 
as lon g as the initial field strength is weaker than B < 10 _5 (n/10 3 cm -3 ) 0,55 G. iMaki fc Susa 
(120071 ) also showed that ionization fraction is high enough for the magnetic field to couple 
with the gas even in the accretion phase after the protostar is formed. In summary, the 
magnetic field is largely dissipated by the Ohmic dissipation before the protostar formation in 
present-day cloud, while the magnetic field can continue to be amplified without dissipation 
in primordial cloud. 

Recent cosmological simulations inferred that a single massive star is form ed without 

: 10 3 cm -3 (Abel et al. 



fragmentat i on in the first collapsed object that is f ormed at n, 



20001 . |2002| ; iBromm et al.l 120021 ; lYoshida et al.l 120061 ) . However, in their simulations, the 
cloud evolution are investigated only for n c < 10 17 cm~ 3 . Since the protostar is formed at 



10 21 cm 3 ( Omukai fc Nishi 19981), the protostar is not yet formed in their simulations. 



On the other hand, 



Machida et al 



(l2007bl ) showed that fragmentation occurs frequently 
after the equation of state becomes hard for n c > 10 17 cm -3 , indicating that the binary or 
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multiple stellar system can be formed also in the early universe. 

The magnetic field in the collapsing cloud is closely related to fragmentation or formation 
of binary or multiple stellar systems. In present-day star formation, the magnetic field 



stron gly suppresses the rotation driven fragmentation (IMachida et al.ll2007d ; iPrice fc Bate 



20071 ). For the primordial cloud, effects of the magnetic field on fragmentation are still 



unknown. 

In this paper, we investigate the evolution of weakly magnetized primordial clouds and 
formation of Population III stars using three-dimensional MHD simulations, and show the 
driving condition of a jet from proto-Population III stars and the fragmentation condition 
in the collapsing primordial cloud cores. The structure of the paper is as follows. The 
framework of our models and the numerical method are given in §2. The numerical results 
are presented in §3. We discuss the fragmentation and jet driving conditions in §4, and 
summarize our results in 55. 



MODEL 



Our initial settings are almost the same as those of IMachida et al.l (l2006dl . l2007dl ). To 
study the cloud evolution, we use the three-dimensional ideal MHD nested grid code. We 
solve the MHD equations including the self-gravity: 



0. 



p— + p(v ■ \/)v 



dt 



dB 

~dt 



VP - — B x(VxB)- pV0, 

47T 

= V x (v x B), 
2 = 4vrG>, 



(1) 

(2) 

(3) 
(4) 



where p, v, P, B, and <fi denote the density, velocity, pressure, magnetic flux density, and 
gravitational potential, re spectively. For gas p ressure, we use the barotropic relation that 
approximates the result of lOmukai et al.l (120051 ) . They consider detailed thermal and chem- 
ical processes in collapsing primordial gas, adopting the a simple one-zone model where the 
core collapses approximately at the free-fall timescale and the size of the core is abou t the 
Jeans length as in the Larson- Penston se lf-similar solut i on (|Penstonlll969l ; lLarsonlll969l ). We 
fit to the th ermal evolution deriy ed by lOmukai et al.l (120051 ) as a function of the density 
(see Fig.l of IMachida et al.ll2007dl ). and the pressure is calculated using the fitted function. 
Therefore, our equation of state is a barotropic. 
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As the initial co nditions of the clouds, we consider the density profile of the critical 



Bonnor-Ebert sphere (jEbertlll955l ; lBonnorlll956l ) with perturbations 



n(r ^ = / PbM / [1 + 6p(r, <p)] for r < R , 

P{ ^ ! I Pbe(Ro) f [1 + 6p(r, ^ for r > R , 1 ) 

where Pbe(^) denotes the density distribution of the critical Bonnor-Ebert sphere. The 
central number density of the Bonnor-Ebert sphere is set at Pbe(O) = 10 3 cm~ 3 . The factor 
/ denotes a density enhancement factor, and the density is increased by a factor / = 1.86 
to promote the collapse. The initial cloud has therefore the central density of no = 1.86 x 
10 3 cm~ 3 . The initial temperature is set at 250 K, and the radius of the Bonnor-Ebert 
sphere is Rq = 6.5 pc, corresponding to the non-dimensional radius of Bonnor-Ebert sphere, 
£max = 6.5. Outside this radius, a uniform gas density of nBE(-Ro) = 130 cm -3 is assumed. 
To promote fragmentation, a small m = 2-mode density perturbation is imposed to the 
spherical cloud such as, 

5p = A^ir/Ro) cos2yj, (6) 

where A v is the amplitude of the perturbation and A v = 0.01 is adopted in all models. The 
total mass contained inside Rq is M c = 1.8 x 10 4 M Q . The density enhancement factor / 
specifies the ratio of the thermal to the gravitational energy a (= -Eth/l-^gravl), where E th 
and -Egrav are the thermal and gravitational energies of the initial cloud). The initial cloud 
has «o = 0.5 when / = 1.86 and A v = 0. 

The initial cloud rotates around the z-axis at a uniform angular velocity of Qq, and it 
has a uniform magnetic field Bq parallel to the z-axis (or rotation axis). The initial model 
is characterized by two non-dimensional parameters: the ratio of the rotational energy to 
the gravitational energy (3 (= -Erot/I-Egravl), and the magnetic energy to the gravitational 
energy 7 (= -E m ag/|-Egrav|) , where E TOt and ^mag are the rotational and magnetic energies. 

We made 36 models changing these two parameters. All the models examined here 
and their parameters (/?o and 70) are summarized in Table [TJ These tabulated values are 
calculated in the cases of spherical clouds, i.e., A v = 0. We examine a large parameter range 
of Pq = 10~ 6 — 0.1, 70 = — 0.2. For convenience, (3q is referred as "initial rotational energy", 
and 70 as "initial magnetic energy." The angular velocity Qq and magnetic field Bq are also 
summarized in Table [TJ In addition, we estimate the mass-to-flux ratio, 

where M and $ denote the mass contained within the cloud radius Rq, and the magnetic 
flux threading the cloud, respectively. There exists a critical value of Mj $ below which a 
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cloud is supported against the g ravity by the magnetic field. For a cloud with an uniform 
density, iMouschovias fc Spitzerl (119761 ) derived a critical mass-to-flux ratio, 



(9. 



3tt V G 



1/2 



where £ = 0.53 for uniform spheres (see also, iMac Low fc Klessenl 12004 ) . For convenience, 
we use the mass-to- flux ratio normalized by the critical value, 



(-) 



t)'S). 



(9) 



which is summarized in Table [TJ 



We adopt the nested grid method (for details, see iMachida et al.l l2005al . l2006al ) to 
obtain high spatial resolution near the center. Each level of a rectangular grid has the same 
number of cells (= 128 x 128 x 64), with the cell width h(l) depending on the grid level 
/. The cell width is halved with every increment of the grid level. The highest level of 
grids changes dynamically: a new finer grid is generated whenever the minimum local Jeans 
length Aj falls below 8h(l max ), where h is the cell width. The maximum level of grids is 
restricted to Z max = 30. Since the d ensity is highes t in th e finest grid, generation of a new 
grid ensures the Jeans condition of iTruelove et al.l (119971 ) with a margin of a safety factor 
2. We begin our calculations with three grid levels (I — 1 — 3). Box size of the initial finest 
grid / = 3 is chosen to be 2R , where Rq is the radius of the critical Bonnor-Ebert sphere. 
The coarsest grid (1 = 1) then has box size of 2 3 Rq. The mirror symmetry with respect to 
z=0 is imposed. A boundary condition is imposed at r = 2 3 Rq, where the magnetic field 



and a mbient gas rotate at an angular velocity of f2 (for deta il, see iMatsumoto &: Tomisaka 
20041 ) . We adopted the hyperbolic V • B cleaning method of iDedner et al.l (120021 ). 



3. RESULTS 

3.1. Typical Jet Model 

Firstly, we show a cloud evolution with a clear jet from a non-fragmentation core. 
Figure [1] shows the cloud evolution before the protostar formation (n c < 10 21 cm -3 ) from 
the initial stage for model 23. Model 23 has parameters of (70, /3 ) = (2 x 10~ 3 , 10~ 4 ), and 
the initial cloud has a magnetic field strength of B = 10~ 6 G, and a angular velocity of 
Qo = 2.3 x 10~ 16 s _1 . Figure [T]a shows the initial spherical cloud (i.e., Bonnor-Ebert Sphere) 
threaded by the uniform magnetic field. Figures (TMH^ show the cloud structure around the 
center of cloud when the central density reaches n c = (6) 4.1 x 10 7 cm -3 , (c) 7.0 x 10 10 cm" 3 , 
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and (d) 5.5 x 10 12 cm -3 , respectively. The density contours projected on the sidewall in 
these figures indicate that the central region becomes oblate as the cloud collapses because 
of magnetic field and rotation, both of which are amplified by the cloud collapse. Figures [TJa- 
[TJi also show that the magnetic field lines are gradually converged toward the center as the 
central density increases. 

Figure [T]e shows the cloud structure at n c = 9.2 x 10 17 cm~ 3 , where a thin disk is 
formed at the center of cloud. In this figure, the disk-like structure is threaded by the 
magnetic field lines that are strongly converged toward the center of the cloud. The density 
contours projected on the sidewall of Figure [T]e show the shock waves above and below 
the disk (see the crowded contour lines on the right sidewall). Near the shock surface, the 
magneti c field lines are bent. Similar structures are seen in the present-day star formation 



process ( IMachida et al.ll2005bl . l2006al ). We discuss the relation between the shock generation 
and the evolution of the magnetic field and rotation in §4.1.31 The cloud structure at 
n c = 3.3 x 10 18 cm~ 3 are shown in Figure [T]f, in which the magnetic field lines are strongly 
converged but slightly twisted. Also in this figure, clear shocks are seen above and below 
the disk-like structure. Figures [Ua-EQf show that the magnetic field lines are hardly twisted 
before the protostar formation (n c < 10 21 cm -3 ), because the collapse timescale is shorter 
than the rotational timescale. 

Figure [2] shows the cloud evolution after the protostar formation (n c > 10 21 cm -3 ). The 
upper and middle panels show the density and velocity distributions on z = and x = 
plane, respectively, and the lower panels show the cloud structures and configurations of the 
magnetic field lines in three dimensions. Figures [2]a-[2]c show the cloud structure 1.4 days after 
the protostar formation. In this model, when the central density reaches n c ~ 4.7 x 10 21 cm -3 , 
a protostar surrounded by the strong shock is formed around the center of the cloud. The 
protostar has a mass of M — 7.6 x 10 -3 M , and a radius of r = 1.1 R Q at its formation epoch. 
The shock front corresponding to protostellar surface is seen in Figure [2]a (r ~ 0.005 AU; 
contours of n ~ 10 21 cm -3 ). The density contours projected on the sidewall in Figures [2]c also 
show the shock surface near the center of the cloud. The magnetic field lines are strongly 
converged to the center of the cloud, but slightly twisted at this epoch. The jet starts to 
be ejected 1.68 days after the protostar formation. Figures [2Je-Eb shows the cloud structure 
2.82 days after the protostar formation. The jet ejected near the protostar is indicated by the 
red contour in Figure [2]f and the red transparent surface in Figure Ejg, which are a boundary 
between the inflow (v r < 0) and jet (v r > 0) regions. The magnetic field lines are twisted 
significantly inside the jet region because of the rapid rotation of the protostar. The jet 
affects the density distribution; the butterfly-like density distribution is shown in Figures [2]f 
and (see contour lines near the protostar projected on the right sidewall). The strong 
shock waves are generated at the upper and lower ends of the jet, reflecting a strong mass 
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ejection. 

Figures [2Fi-T2]/ show the cloud structure 3.48 days after the protostar formation, the 
last stage of this calculation. Until the last stage, the jet continue to extend as is seen 
in Figures [2]f and [2]j. In Figure [2pi, two nested shocks are observed near the center of 
the cloud (r ~ 0.005 AU and r ~ 0.01 AU). The outer shock is generated by the mass 
ejection (i.e., jet) near the protostar, while the inner one corresponds to the surface of the 
protostar. In Figure [2]/, the isodensity surface exhibits a shallow cone-shape on the disk, 
which corresponds to the outer shock in Figure [2Fi. The magnetic field lines inside the jet 
region are more strongly twisted with times as shown in Figures Wq and [2]/. At the end 
of the calculation, the jet has a maximum speed of t> max = 66.3 km, and it extends up to 
0.05 AU. The mass of protostar reaches M ps = 8.5 x 1O _3 M , while mass of the outflowing 
gas is M out = 1.1 x 1O~ 3 M . About 10% of the accreting matter is therefore ejected from 
the center of the cloud via the protostellar jet. The mass ejection rate is estimated as 
M = 5.9 x lO^Moyr- 1 . 



3.2. Typical Fragmentation Model 

Next, we describe the cloud evolution exhibiting fragmentation. Figure [3] shows the 
cloud evolution for model 21 with (70, flo) — (2-5 x 10" 7 , 10~ 4 ), which are corresponding to 
a magnetic field strength of Bq = 10~ 8 G and an angular velocity of Qq = 2.3 x 10 _16 s~ 1 
at the initial stage. Model 21 has the same rotational energy as the model 23, the previous 
model, but 10~ 4 times smaller magnetic energy than model 23. Figures [3]a and[3]6 show the 
cloud structure just before the protostar formation (n c = 1.0 x 10 17 cm -3 ). We define the 
protostar formation epoch as the stage of n c = 10 21 cm" 3 . Similar to model 23, a thin disk is 
formed at the center of the cloud before the protostar formation. The disk keeps an almost 
axisymmetric structure in the plane perpendicular to the rotation axis. Although we added 
1% of the non-axisymmetric density perturbation to the initial state, the non-axisymmetric 
structure hardly evolves before the protostar formation. Figure [3f> also shows an hourglass 
configuration of the magnetic field lines, similar to model 21. 

Figures [3]c and[3]<i show the cloud structure 2.4 days after the protostar formation. The 
central region is deformed from a disk structure into a ring at n c ~ 8 x 10 18 cm" 3 , and th e 



ring fragments. This is consistent with the fragmentation criterion of lMachida et al.l (j2007dl ); 
a non- magnetized cloud with rotational energy of /3q > 10" 5 undergoes fragmentation. The 
cloud here has a rotational energy of /3q = 10 -4 and the very small magnetic energy of 70 = 
2.5 x 10~ 7 at the initial stage. Fragments are located at (2, y) ~ (±0.035 AU, =f0.055AU) 
with long spiral tails that are the remnant of the ring (Fig. [3]c) . The magnetic field lines are 
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distributed strongly along the spiral tails (Fig. [3]rf). 

Figures [3]e and [3]f show the cloud structure 8.8 days after the protostar formation. The 
fragments exhibit a round shape in Figure [3]e. The central density of the fragments reaches 
n c ~ 1.8 x 10 22 cm -3 at this epoch. Fragments are surrounded by their respective disk with 
spiral arms, and the disks grow with time, comparing Figures |3]c and [3]e. At the end of 
the calculation, each fragment (i.e., each protostar) has a mass of ~ 1.3 x 1O~ 2 M and a 
radius of ~ 1.4 R Q . The separation between fragments is ~ 0.89 AU. As shown in Figure [3]f, 
the magnetic field lines are distributed along the spiral tails, indicating that the vertical 
component of the magnetic field is weaker than the other component (i.e., B z <C B x ,B y ). 
This is responsible to the weak magnetic field around the protostars, and the magnetic field 
lines easily follow the orbital motion of the fragments. 

This fragmentation model does not exhibit a jet formation by the end of the calculation 
in spite of the same initial rotation speed as in the previous jet model. This is also responsible 
to the weak magnetic field. Note that the calculation halts 8.8 10 days after the stellar core 
formation, and we does not reject a possibility that a jet appears in later stage. 

One may image that this model undergoes the rotation driven fragmentation because of a 
weak magnetic braking, and the previous model is kept away from the fragmentation because 
of the strong magnetic braking. However, both the models show the same amplification rate 
for the magnetic field and the angular velocity during the collapse. The saturation levels of 
the amplification for the magnetic field and angular velocity depend on the initial conditions. 
The two models track two different evolutional sequences according to their respective 70/A). 
This point is discussed in §4.21 and $51 

3.3. Fragmentation and Magnetic/Rotational Energy 

Figure H] shows the final stages in the plane of initial magnetic (70, x-axis) and rotational 
energies (/?o, y-axis) for every models. The density distribution in the z = plane is shown 
in each panel. Models are classified into four types: fragmentation models (models in the 
blue background), non- fragmentation models (models in the red background), merger models 
in which the fragments merge after fragmentation (models in the green background), and 
no collapsing model (model in the gray background). In fragmentation models, calculations 
were stopped when the Jeans condition was violated in a grid except for the finest grid. Such 
a case occurs either when fragments escape from the finest grid or when the gas far from 
the center becomes denser than that in the finest grid. In non-fragment models, the Jeans 
condition was not violated because the densest gas is located at the finest grid, and the 
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calculations was stopped at some high density. After the protostar formation, the Alfven 
velocity becomes extremely large due to the strong magnetic field at the central region, 
and the timestep becomes extremely short. We stopped calculation after we checked that 
fragmentation was not likely to occur around the protostar. The fragmentation reproduced 
here is therefore restricted to the cases where fragmentation occurs near the central region 
and fragmentation occurs just after the protostar formation. In merger models [(70, (3q) = 
(2 x 10~ 3 , 10~ 2 ), (2 x 10~ 3 , 10~ 3 )], fragments merge to form a single core at the center of the 
cloud after fragmentation. The merged core did not undergo fragmentation again although 
we calculated merger models for sufficiently long term. In the no collapsing model [(70, (3q) 
= (0.2, 0.1)], the cloud oscillates without collapse because this model has the large magnetic 
and rotational energies. 

In Figure HI fragmentation models are distributed in the upper left region, indicating 
that a cloud tends to fragment when the initial cloud has a weaker magnetic field and faster 
rotation. Namely, the magnetic field suppresses fragmentation, while the rotation promotes 
fragmentation. This effect of the magnetic field and rotation on fragmentation is similar to 
that in present-day star formation (Hoski ng fc Whitworthl 12004 ; iMachida et al.ll2004l . l2005bl . 
2007cl : IPrice fc Batel 120071 : Irlennebelle fc Tevssier! 1 20071 ). Firstly, to investigate the effect of 
rotation on fragmentation, we focus on models displayed in the fourth column in Figure H] 
(models 4, 10, 16, 22, 28, and 34). These models have the same magnetic energies of 70 = 
2x 10~ 5 and different rotational energies of (3q = 10~ 6 — 0.1 at the initial stage. Fragmentation 
occurs in models with (3 Q > 10 -4 (models 4, 10, 16, 22), while fragmentation does not occur 
in models with (3 < 10 -4 (models 28, 34). In addition, models with larger rotational 
energies undergo fragmentation in the earlier evolution phases with wider separation. For 
example, model 4 (/?o = 0.1) fragments at n c = 9.5 x 10 17 cm~ 3 with the separation of 
R sep ~ 4.3 AU, while model 22 (0 n = 1 0~ 4 ) does at n c = 1.6 x 10 19 cm -3 with the separation of 
R sep ~ 0.2 AU. IMachida et al.l (l2007dl ) also showed that a cloud with faster rotation exhibits 
wider separation. On the other hand, paying attention to the non-fragmentation models, 
model 34 exhibits a compact spherical protostar with a disk, while model 28 shows bar and 
ring structures. Model 34 did not show any sign of fragmentation, and the central core 
continues to grow. Model 28 has marginal parameters for fragmentation since it is located 
at the boundary between fragmentation and non-fragmentation models. The protostar fails 
to fragment in spite of its bar shape. The spiral arms evolve into a ring, which does not also 
fragment. 

We now focus on models displayed in the fourth low in Figure H] (models 19-24), to 
investigate the effect of magnetic field on fragmentation. These models have different mag- 
netic energies of 70 = — 0.2 but the same rotational energies of (3q = 10 -4 at the initial 
stage. In these models, weakly magnetized clouds undergo fragmentation (70 < 2 x 10~ 5 , 
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models 19-22), while strongly magnetized clouds do not (70 > 2 x 10~ 5 , models 23 and 24). 
When clouds have the same rotational energies at the initial stage, model with weaker mag- 
netic field undergo fragmentation at the earlier evolutionary phases with wider separation. 
For example, the model with weaker magnetic field 70 = 2 x 10 -9 (model 28) fragments at 
n c = 7.8 x 10 18 cm~ 3 , while the model with stronger magnetic field 70 = 2 x 10 -5 (model 
22) does at n c = 1.6 x 10 19 cm" 3 . Moreover, the models with much larger magnetic energies 
(models 23 and 24) produce a single compact core at the center of the cloud. 

Fragmentation and non-fragmentation models are clearly separated in Figure HI In 
addition, merger models are located at the boundary between fragmentation and non- 
fragmentation models. In fragmentation models, models with weaker magnetic field and 
faster rotation tend to have wider separation. In non-fragmentation models, clouds have 
more compact cores in models with stronger magnetic fields and slower rotations. These fea- 
tures indicate clearly that the rotation promotes fragmentation, and the magnetic field sup- 
presses fragmentation in primordial star formation, as well as the present-day star formation. 
The models presented here have a small amplitude of the non-axisymmetric density pertur- 
bation (A<p = 0.01) at the initial stage. If larger amplitude was set at the initial stage, the 
boundary between the fragmentation and non-fragmentation models might change slightly. 



Machida et al.l (j2007d ) discussed the fragmentation condition for a present-day magnetized 



cloud, and indicated that the fragmentation condition does not dep end qualitativel y on the 
initial amplitude of the non-axisymmetric perturbation. Moreover, iMachida et al.l (l2007dl ) 
showed that fragmentation condition for a primordial non-magnetized cloud depend slightly 
on the initial amplitude of the non-axisymmetric perturbation. We therefore expect that 
the initial amplitude of the non-axisymmetric perturbation hardly affect the fragmentation 
condition even in the primordial magnetized collapsing cloud. 



3.4. Jet and Magnetic/Rotational Energy 

Figure shows final states in the plane of the initial magnetic (70, x-axis) and ro- 
tational energies (/3q, y-axis), contrasting the jets at almost the same evolutionary stage of 
n c ~ io 22-23 cm -3 . Each panel shows the density and velocity distribution in the y = plane, 
and a thick red contour denotes the boundary between inflow (v r < 0) and outflow (v r > 0) 
regions. The gas flows out of the central region inside the red contour. The panels without 
thick red contours indicate the models where a jet does not appear. As listed in Table [TJ jets 
appear in 13 of 36 models. We call these 13 models 'jet models', and other models 'non-jet 
models.' In all the jet models, a jet appears only after a protostar is formed. After the pro- 
tostar formation, the rotation timescale becomes shorter than the collapse timescale because 
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of the accretion of the angular momentum and hardness of the equation of stage at the cloud 
center, and the strong centrifugal force drives the jet via disk wind mechanism ( ITomisaka 



20021 ; iBaneriee fc Pudritzll2006l : Machida et al.ll2007bl ; iHennebelle fc Fromandl2007n . The jet 



disturbs the density distribution near the protostar. As shown in Figure the jet supplies 
gas above and below the disk, and the density becomes higher there for the jet models, while 
these region remains less dense for the non-jet models. 

The jet models are distributed in the lower right region in Figure 0, indicating that a 
jet is driven in strongly magnetized but slowly rotating cloud. It should be also noted that 
the non-jet models coincide with the fragmentation models except for models 31 and 32. In 
other words, almost all the clouds experience either jet formation or fragmentation. In the 
fragmentation models, the angular momentum of the parent cloud is distributed to both the 
orbital and spin angular momenta of the fragments, and fragmentation therefore reduces the 
spin angular momentum of the protostar. This deficiency of the spin angular momentum of 
the fragments suppresses jet formation. 

Firstly, to investigate the relation between jet and magnetic field strength, we focus on 
models in the fourth row (models 27-30). In model 27, which has a weakest magnetic field 
B = 1CT 8 G (or 70 = 2 x 1CT 7 ) at the initial stage, no jet appears because of fragmentation. 
In model 28, the protostar manages to drive the jet, while the initial cloud has a considerably 
small magnetic energy (/?o = 2 x 10~ 5 ). When the cloud has a weak magnetic field, a jet 
is considered to be driven by the magnetic pressure gradient force, not by the disk wind 
mechanism. The j et structure in mode l 28 is simil ar to the magnetic bubble as rep r oduced by 
Tomisakal (120021 ) , iBanerjee fc Pudritzl (120061 ) and IHennebelle fc Fromang! (120071 ) . ITomisaka 
( 20021 ) shows that when the magnetic field around the central object is extremely weak 
(/? p ^> 1, where (3 P is the plasma beta), the magnetic field is amplified by the spin of the 
central object, and the magnetic pressure drives a jet. 

A prominent jet is reproduced in model 29, in which the jet is considered to be driven 
mainly by the disk wind mechanism. The jet in model 30 is weaker than that in model 29 
in spite of the stronger initial field strength; the jet in model 30 has slower average speed 
and smaller mass ejection rate than those in model 29. The weakness of the jet in model 
30 is attributed to the strong magnetic braking, which slows down the spin of the cloud 
before the protostar formation in model 30. In fact, the protostar has an angular velocity 
of f2 c = 1.0 x 10 _5 s _1 in model 29 while fi c = 1.1 x 10~ 6 s _1 in model 30 at the protostar 
formation epoch. The protostar therefore rotates ~ 10 times faster in model 29 than that 
in model 30. In addition, the growth rate of the angular velocity in the collapsing cloud 
is smaller in strongly magnetized cloud (see §4. 1.11) . In models 27-30, the jet in model 29 
has highest speed and largest jet momentum flux at the end of the calculation. As a result, 
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moderate strength of the magnetic field is necessary to drive strong jet from the protostar. 

Next, to investigate dependence on rotation, we focus on the models in the third column 
in Figure (models 11, 17, 23, 29, and 35). In these models, the most prominent jet appears 
in model 23, while weak jets appear both in slowly (model 35) and rapidly rotating clouds 
(model 11). The weakness of the jet in model 35 is due to deficiency of rotation to drive a 
strong jet, and that in model 11 is due to a low amplification of the magnetic field during 
the collapse; amplification of the magnetic field is related to the rotation, and spin-up is also 
related to the magnetic field strength as discussed in §4.1.11 Therefore, a moderate rotation 
speed as well as a moderate magnetic field strength is necessary to drive a strong jet. In 
summary, a jet appears in models distributed in the lower right region (i.e., models with 
stronger magnetic field and slower rotation), and models showing strong jets are limited in 
the range of 1(T 5 < /3 < 1(T 3 and 1(T 3 < 7 < 10 _1 . 

Since we calculated the evolution of the jets only for short duration (~ a few 10 days) for 
numerical limitation, we cannot estimate the momentum of a jet, mass ejection speed and 
final speed of jet in long-term evolution. In addition, strength of a jet may largely change 
in further long-term evolution. We do not discuss the further evolution of a jet (c.f., final 
jet speed, mass ejection rate, etc..) , because our purpose in this paper is to investigate the 
condition for driving jet in primordial collapsing clouds. 



3.5. Configurations of Jets and Magnetic Field Lines 

Figure shows a configuration of magnetic field lines (black-white streamlines) and a 
structure of a jet (transparent red iso-velocity surface) for each models. 

In the jet models, the magnetic field lines are twisted inside the jet region. Almost 
axisymmetric jets are seen in models 18, 23, 24, 28, 29, 34, and 35. These jets have the 
hourglass-like configurations of the magnetic field lines, where the poloidal component is 
more dominant than the toroidal component as shown in Figure [T|f and [2]c. These con- 
figur ations of the magnetic fie ld lines can easily drive a jet by the disk wind mechanism 



(e.g., iBlandford fc Payndll982l ). Non-axisymmetric jets are seen in models 17, 30, and 36. 
They are caused by the non-axisymmetric density distribution at the protostar formation 
epoch. Such a non-axisymmetric density distribution is also reflected by non-axisymmetric 
circumstellar disks, represented by a red iso-density surface in Figure El For example, model 
30 exhibits the non-axisymmetric jet (transparent red iso-velocity surface) driven from the 
bar-like density distribution (red iso-density velocity surface) at the root of the jet. 



In the fragmentation models, the configurations of the magnetic field lines are disturbed, 
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and the toroidal field tends to be more dominant than the poloidal filed as shown in models 
16 and 22 in Figure [6] and also in Figure [3]f. These configurations of the magnetic field 
lines indicate weakness of the field strength, which is insufficient to drive a jet. The weak 
poloidal field is attributed to the oblique shocks on the surface of the disk envelope, and the 
disturbance of the field lines are attributed to the orbital motion of the fragments and the 
spin of the fragments. The spin of the fragment winds the magnetic fields around its rotation 
axis, and the field strength is amplified. In the further stages, the magnetic pressure may 
drive a jet. 



4. DISCUSSION 
4.1. Magnetic Flux-Spin Relation 

4-1.1. Evolution Track for Magnetic Field and Angular Velocity 

Figure [7] shows the evolution track of the central magnetic field strength (x-axis) and 
central angular velocity (y-axis) from the initial state for some models. The x-axis indicates 
B c /(8tcP c ) 1 ^ 2 , which corresponds to the pressure ratio of (P mag / P c ) 1 ^ 2 at the center, where 
P mag and P c denote the central magnetic pressure and the central thermal pressure. This 

— 1/2 

variable also corresponds to the inverse square root of the plasma beta, j3 v . Hereafter this 
variable is called the normalized magnetic field. The y-axis indicates Vt c / (AnGp^) 1 ! 2 , which 
corresponds to a ratio of the time scales, t s /t rot , at the center, where t s and t rot denote 
the free-fall and rotation time scales. This variable also coincide with the energy ratio of 

1 /2 

(_E ro t/|-Eg r av|) at the center if assuming a uniform sphere. Hereafter this variable is called 
the normalized angular velocity. The asterisks at the ends of the loci denote the initial values 
of normalized magnetic fields and normalize angular velocities. The diamonds on the loci 
denote the stages of n c = 10 5 , 10 7 , • • •, 10 19 cm -3 , and these values are labeled for the loci of 
models 2, 9, 21, 23, 33, and 35. The arrows on the loci indicate the directions of the cloud 
evolutions. 

The normalized magnetic fields and angular velocities increase as the central density in- 
creases except for stages of very high density in all the models. For comparison, a relationship 
of 

(4vrG Pc )V2 K (8ttP c )V2- (10) 

is plotted with the solid line at the lower right corner of Figure [TJ and all the model evolves in 
the direction almost parallel to this line. This implies that both of the normalized magnetic 
fields and the normalized angular velocities are amplified during the early phase of the 
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collap se at the almost same growth rate. As indicated by iMachida et al.l (l2005al . [2006a 



2007al ). the central magnetic field and angular velocity increase as, 

b c , n c oc p 2 J 3 (ii) 

during the nearly spherical collapse phase that corresponds to Figure [T]a-[TK The relationship 
of the angular velocity in equation ( flQj) is rewritten as, 

(isfe* K (12) 

and the relationship of the magnetic field is rewritten as, 

— ^ - (13) 



57T 



The power index in the right hand side of f[T3"j) is slightly lower than that in equation (fT2"l) 
although they coincide in the isothermal case. The difference between the power indexes 
is attributed to the equation of state in the primordial gas, P = p 1 . The polytropic index 
7 is assumed as a function o f density, and i t is ap proximated to 7 ~ 1.1 for 10 4 cm -3 < 
n < 10 17 cm~ 3 (see Fig.l of IMachida et al.l l2007al ). Adopting the relation of P c oc p^ 1 , 



the equation (1131) is obtained. The slightly larger power index in the right hand side of 
equation ffl3|) is reflected in the slope of loci slightly steeper than that of equation fflUj) . The 
relationships of equations ffl2l and ffl3|) are examined in detail in §4.1 .21 

The slope of loci also indicates that the magnetic braking is ineffective during the 
spherical collapse phase. If magnetic braking is effective (i.e., it slows down the spin of a 
cloud), increase in the angular velocity was less than that indicated by equation ([121) . and 
the slope of the loci would be shallower. We confirmed that the magnetic braking become 
prominent only after the protostar formation. In other words, the magnetic field is not 
amplified enough to affects the evolution of the cloud during the spherical collapse of the 
early phase. In the later phase, the cloud center exhibits transition from nearly spherical 
collapse to the disk-like collapse because of rapid rotation and/or strong magnetic field, which 
are amplified in the nearly spherical collapse phase, and the evolutionary track stagnated in 
the hatched band of Figure UJ The disk-like collapse phase is discussed in §4.1.31 



4-1.2. Magnetic Field and Angular Velocity in the Spherical Collapse Phase 

To verify the relation of equations ffl2|) and ffl3|) . the evolutions of the central angular 
velocities and magnetic fields are plotted against the central density in Figures [8] and [91 
respectively. 



-16- 



Figure [8] shows the evolution of the normalized angular velocity in models 3, 9, 15, 21, 
27, and 33, which are located in the third column of Figure H] and have the common weak 
magnetic field B = 10 -7 G and different angular velocities f2 = 2.3 x 10~ 17 — 7.2 x 10~ 15 s _1 
at the initial stage (see Table [1]). In Figure El the normalized angular velocities in all 
models except for model 3 evolve with oc p 1//g from the initial stage, indicating the power 
law amplification of fi c oc p 2 J Z and spherically collapse of the clouds. On the other hand, the 
cloud in model 3 cannot collapse spherically because the cloud has a large rotational energy 
at the initial stage. 

Figure [9] shows the evolution of the magnetic field strength normalized by the square root 
of the density (Fig. [9] upper panel) and the normalized magnetic field (Fig. [9] lower panel) in 
models 32-36. Models 32-36 are located in the sixth row in Figure HI and have the common 
small angular velocity Qq — 2.3 x 10 _17 s _1 and different strengths of the magnetic field 
Bq = 10~ 9 — 10~ 6 G at the initial stage (see Tabled]). The upper panel indicates the power 
law amplification of B c oc p 2 ^ 3 (equation [IT] ) except in the high-density phase and therefore 
the spherical collapse of the clouds for all models. Figure [9] lower panel also shows the 
evolution of the magnetic field in terms of the j3 p (= B c j \/87rP c ). All the curves indicate 
the power law amplification of oc p ' 12 , corresponding to equation (113j) . The curves in 
the lower panel undulate more than those in the upper panel because of the undulation in 
the pressure included in the normalization. However, the curves in the lower panel converges 

— 1/9 

well to j3 p ~ 0.2 at the high density phase compared with those in the upper panel. As a 
result, when the cloud collapses spherically (i.e., early evolution phase), the evolution of the 
angular velocity and magnetic field are well expressed by equation (|lip . and the normalized 
variables are also expressed by equations (fl2"l) and (fl3"|) . 



4-1.3. Convergence of the Magnetic Field and Angular Velocity 

As the cloud collapses, the normalized angular velocities converge to il c /(4:7iGp c ) 1 ^ 2 ~ 
0.25 in Figure [HI and the normalized magnetic fields converge to B c /(8irP c )V 2 ~ 0.2 with 
initially strong magnetic fields (models 34, 35, and 36) in Figure [9] lower panel. Figures [8] 
and [9] present only limited models. As shown in Figure [7] where more models are plotted, 
the normalized magnetic field and the normalized angular velocity converge to the so-called 
flux-spin relation, 

O 2 R 2 

^ + ^ = 1. (14) 

(0.25) 2 47rGp c (0.2) 2 8vrP c 1 ; 

This relation is represented by the gray band in Figure [71 If an initial stage is located in 
lower and left-side of the gray band in Figure [7J a points indicated by the two normalized 
variables moves in the upper right direction according to equations f[T2l and ([TBI in the 
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nearly spherical collapse phase. When the point reaches the gray band of equation ()T4j) . the 
cloud changes the geometry of collapse from a sphere to a disk. After the transition, the 
point remains around the gray band with small oscillation. 

The numerical factors in equation (Tl4"l) are determined empirically. The horizontal gray 
band corresponds to the convergence value in Figure E [fi c /(4vrG'p c ) 1 / 2 ~ 0.25], and the verti- 
cal gray band corresponds to the convergence value in Figure [8] lo wer panel \B r J {S-kPc) 1 ^ 2 ~ 



0.2]. The convergence value of the ang ular velocity was f ound by iMatsumoto et all (119971 ). 



and that of magnetic field was found by lNakamura et all (119991 ). for isothermal clouds mod- 



eling the present-day star form ation. While they considered either rotation or magnetic field, 



Machida et all (j2005al . l2006al ) considered both the rotation and magnetic field, and obtained 



the convergence values for rotation and magnetic field (flux-spin relati on) for the present-day 



star fo rmation. Equation (|14[) is an extension of flux-spin relation of iMachida et al.1 ((2005a 



2006al ) to the primordial star formation. Figure [7] indicates that the magnetic field strengths 
and angular velocities converge to the flux-spin relation represented by equation (THl) even 
in the primordial magnetized clouds. 

The length of the locus before the convergence to the spin-flux relation depends on 
the initial magnetic field strength and the initial angular velocity. Figures [8] and [9] show 
that models having stronger magnetic field and/or more rapid rotation at the initial state 
converge to the spin- flux relation in the earlier evolutionary phase (or lower density phase). 
In other words, a model initially located near the gray band converges earlier. For example, 
in Figure [3, the evolutionary track of model 9 reaches gray band at n c ~ 10 9 cm~ 3 , while 
that of model 21 reaches the gray band at n c ~ 10 13 cm~ 3 . This is because the initial 
normalized angular velocity in model 9 is stronger than that in model 21 (models 9 and 21 
have the same strengths of the magnetic field at the initial). In order to confirm the relation 
of equations (1T21 and (TIBl in Figure [3, we focus on the evolutionary track for model 33, 
which has fi c /(47rGp c ) 1 / 2 = 5.5 x 10~ 4 and B c / (8irc 2 s p c ) = 1.4 x 10~ 4 at the initial stage. 
The evolutionary track in model 33 reaches gray band at n ~ 10 19 cm -3 , that is 10 16 times 
higher than the initial density. According to equations (TT2|) and ([TBI , the normalized angular 
velocity and magnetic field should be amplified by factors of 464 and 83.2 respectively. On 
the other hand, the model 33 exhibits fi c /(47rGp c ) 1 / 2 ~ 0.2 and BJ (Mc 2 s p c ) = 1 x 10~ 2 
when the evolutionary track reaches gray band, and the amplification factors are therefore 
364 and 71.4, which are consistent with that derived from equations (fl2l and ([TBI . 

When an evolutionary track converges to the flux-spin relations, the central cloud de- 
forms from a sphere to a disk because the central angular velocity and/or magnetic field 
have been amplified during the nearly spherical collapse phase. For example, model 23, 
shown in Figure (U produces a thin disk at n c ~ 10 13 cm -3 , at which the evolutionary track 
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converges to the gray band in Figure [7J After evolutionary track converges to the flux-spin 
relation, it oscillates around the flux-spin relation as shown in Figures [7J We also noticed 
that the strong accretion shocks are generated on the surfaces of the disk when the model 
converges to the flux-spin relation. The shock generation is also related to the oscillation of 
the evolutionary track after the convergence. When the evolutionary track oscillates around 
the flux-spin relation, a new shock is generated inside the collapsing disk intermittently. 
The intermittent shock generation results in a nested structure of the shocks. Such nested 
shocks are firs t repr oduced in the simulation of collapse of a rotating isothermal disk by 
Norman et al.l (Il980l ). The convergence and oscillation is also seen in Figures [8] and [9j 



The convergence indicates a power law amplification of angular velocity and magnetic 
field during the disk-like collapse phase. In the case of homologous collapse of an isothermal 

1/2 

disk with a thin disk ap proximation, this power law amplification is expressed as f2 c , B c oc p c 
( Machida et al. 2005a ). For the polytropic gas, the growth rate is modified as Q c , B c oc p 7 / 2 . 



The normalized va riables are also modifi ed as fi c /(47rGp c ) 1 / 2 oc pi 7 and _B c /(87rP c ) 1/2 



constant (see also iMachida et al.ll2007dl ). In our case of 7 m 1.1, the normalized angular 
velocity would increase slightly in proportion to p^' 05 , an d such a small increase is hardly 
seen in Figures [7J and M because of the significant amplitude of the oscillation. 



4.2. Fragmentation/ Jet Condition 

All the models examined in this paper exhibit convergence to the flux-spin relation of 
equation ( fT4l) as shown in Figure UJ A convergence point within the flux-spin relation de- 
pends on the initial angular velocity and the initial magnetic field, and the convergence point 
characterizes the fate of a cloud. A model with initially fast rotation and weak magnetic 
field starts its evolutionary track at a point in the upper left region, and the evolutionary 
track converges to the horizontal gray band. Such a model is called "the rotation-dominated 
model." On the other hand, a model with initial slow rotation and strong magnetic field 
starts its evolutionary track at a point in the lower right region of Figure [7J and the evo- 
lutionary track converges to the vertical gray band. Such a model is called "magnetically- 
dominated model." We classified all the models into the rotation (Q) and magnetic (B) 
dominated models and summarized it in the last column of Table [TJ 

Only the rotation-dominated models exhibit fragmentation because the rotation is am- 
plified enough to promote fragmenta tion and magnetic braking is insufficient to prevent 



fragmentation. IMachida et al.l (l2007dl ) examined fragmentation of non-magnetized primor- 



dial clouds, and showed that fragmentation occurs in the case where the normalized angular 
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velocity at the center converges to 

Q, 



c 



0.25 (15) 



before the protostar formation. A magnetically-dominated model never converges to equa- 
tion (fT5"I) because the evolutionary track never reaches the horizontal band in Figure [7J A 
rotation-dominated model is satisfied with equation (TI5T) when the initial normalized angular 
velocity is 

Q ° >0.25^V /6 , (16) 



y/AnGpo ~ VPcr/ 

where p cr ~ 10 17 cm -3 denotes an upper-bound of density where the equation state is ap- 
proximated by the polytrope with 7 rs 1.1, and po — 1-86 x 10 3 cm -3 denote the initial 
central density. We used the relationship of fi c oc p 2 J z when equation (Fl6|) is derived. Ex- 
amining the initial conditions of the rotation-dominated models, all the rotation-dominated 
models except for models 31, 32, and 33 satisfies equation (Tip]) , and these models exhibit 
fragmentation, consistently with the fragmentation condition. The fragmentation condition 
of equation ffl5|) is therefore valid even for magnetized primordial clouds. 



The fragmentation condition can be rewritten in terms of the rotational energy. As- 
suming that the center of the cloud undergoes homologous collapse with rigid rotation, the 
parameter j3{= E rot /\E zrax \) is expressed as j3 = fl 2 / '(4irG p) . Equation (fT5l) therefore indi- 
cates that fragmentation occurs when the rotational energy reaches ~ 6% of the gravitational 
energy owing to spinning up of the cloud center. 

All the magnetically-dominated models presented here show the jet formation. The 
condition for jet formation coincides with the flux-spin relation of the magnetic-dominated 
models, 

- ~ 0.2. 17 

When the central magnetic field is amplified up to the level of equation (TIT]) before the 
protostar formation, a jet appears just after a protostar is formed. This condition is expressed 
as, 

d / \ °- 12 

B ° > 0.2 (*] , (18) 



using the power law amplification of the magnetic field during the spherical collapse (eq. 
[TB"]). Examining the initial conditions of the magnetically-dominated models, they are 
satisfied with equation (Tl8|) except for model 34, although the evolutionary track of model 
34 converges to the point satisfying equation fflTj) . This exceptional feature of model 34 is 
attributed to the slightly steeper slope in the lower panel of Figure[8]than that of the reference 
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line, (3 P oc n° c 12 . Equation f fTTl) can be rewritten as (3 V ~ 25 since B c / (RitP) 1 / 2 = (3 P 1 ^ 2 . 
This indicate that a jet is driven by the protostar when the magnetic field is amplified and 
the plasma beta reaches f3 ~ 25 before the protostar formation. The equation (|T7|) cannot 
be realized in rapidly rotating clouds, because the evolutionary track converges to horizontal 
gray band in Figure [7] and never reach the vertical band. 

Although the magnetically and rotation dominated models are separated clearly, some 
models are located on the border between the categories present mixed features. Model 11 
and 17 are located near the border between the rotational and magnetic dominated models, 
and exhibit both fragmentation and jet. In these models, fragments merge to form a single 
core after fragmentation, while other fragmentation model does not show merger until the 
end of the calculation. After the merger, the jet begins to be driven from the merged core. 

In the fragmentation models, we do not observe a jet driven from either fragment. This 
indicates that fragmentation occurs in rotation dominated models where the magnetic field 
is too weak to drive the jet. However, it may be possible that the magnetic field amplified by 
the spin of the protostars drive the jet in the further stage. We conclude that the condition 
of equation ffl7|) is valid for the jet which is formed within a few times 10 days after the 
protostar formation. 



4.3. Magnetic Fields and Rotation Periods of Proto-Population III Stars 

In our calculation, both the normalized angular velocity and normalized magnetic field 
are bounded by values expressed by equations ffl~5l) and f|T7|) (i.e., gray band in Fig. [7]). 
Therefore, given a central density, the maximum angular velocity and magnetic field strength 
are expressed as, 

Onax < 0.25v/4vrGp c , (19) 

and 

5 max < 0.2 y/toP e . (20) 

When the protostar is formed at n c = 10 21 cm" 3 , which is 10 4 times larger than n cr , the 
protostar has a maximum angular velocity of 

n ps = 1.1 x 10~ 5 s~\ (21) 

which corresponds to the rotation period of P > 6.6 days. The maximum magnetic field 
strength is given by, 

B ps = 2.6 x 10 4 G. (22) 
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at the protostar formation epoch (n = 10 21 cm~ 3 ), where the thermal pressure adopted in 
our calculation is used. We measure the magnetic field strengths (B ps ) and angular velocities 
(f2 ps ) at the protostar formation epoch for every models, and the measured values are listed in 
7th and 8th column of Table HJ For the models with parenthesis, we cannot follow the cloud 
evolution by the protostar formation epoch because fragmentation occurs in relatively early 
phases. For these models, the magnetic field and angular velocity just before fragmentation 
are listed. As listed in Table (H the protostars at their formation epoch have the angular 
velocities in the range of 3.3 x 10 _7 s^ 1 < Q ps < 1.1 x 10 _5 s _1 , which corresponds to the 
rotation periods of 6.6 days < P < 220 days. These angular velocities are well bounded by 
the value of equation (T2~Tj) . The magnetic fields are ranged in 0.01 G < B ps < 3.1 x 10 4 G, 
which are also bounded by the value of equation (1221) . 

The observati ons indicate that the pres en t-day protostar have magn etic fields of ~ 1 kG 
at the maximum Jjohns-Krull et al.lll999aH bl. 1 200 ll ; iBouvier et al.l 120061 ) . According to our 
simulations, Population III protostars can possess 10 times stronger magnetic field than 
protostars at present day. The smallness of the magnetic field in the present day protostars 



i s attri buted to the magnetic dissipation during the protostellar collapse. iMachida et al. 
(j2007al ) studied collapse of magnetized clouds for the present-day star formation, and showed 
that the magnetic fields of the protostar is ~ kG at the maximum because the magnetic field 
is largely dissipate d by the Ohmic dis sipation in the late phase of collapse, 10 11 < n < 



10 cm (see also. lNakano et al.ll2002l ). On the other hand, the dissipation of the magnetic 



field is not s upposed to be effective in primordial collapsing clouds as shown in lMaki &: Susa 
( 120041 . 120071 ) . Population III stars therefore possess stronger magnetic field than that in 
present day if they are formed via the process of the magnetically-dominated models. 

The magnetic fields and angular velocities of the protostar listed in Table H] are the 
values at the moment of the protostar formation, in which the mass of protostars is only 
~ 10~ 3 — 10 _2 M Q . Since stars acquire a large fraction of mass in a subsequent accretion 
phase, the magnetic field strength and angular velocity may be changed in further evolution. 
The angular momentum is transferred by magnetic interaction between the protostar and 
circumstellar disk. The magnetic field can be amplified by convection inside the protostar. 
However, since the purpose of this paper is to investigate the magnetic effect in collapsing 
primordial clouds, we do not discuss subsequent evolution of the magnetic field and angular 
velocity. To determine the magnetic field and angular velocity of Population III stars, further 
long-term calculations is necessary including a model of stellar evolution. 
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5. SUMMARY 

In this paper, we calculated cloud evolution from the stage of n c = 10 3 cm -3 until 
the protostar is formed (~ 10 22 cm -3 ) for 36 models, parameterizing the initial magnetic 
field strength and the initial rotation, to investigate effects of magnetic fields in collapsing 
primordial clouds. In figure [TU1 the fates of all the clouds are plotted in the plane of the 
parameters 70 and (3q, where the circles, squares, triangles, and crosses mean the models 
showing fragmentation, jet, both fragmentation and jet, and neither fragmentation nor jet, 
respectively. In model located in the upper right corner (model with diamond), the cloud 
oscillates around the initial state without collapse, because strong magnetic field and rapid 
rotation suppress the collapse of the cloud. The solid line in Figure [TU1 represents /3 = 70, 
indicating that the magnetic energy is more dominant than the rotational energy in the 
model distributed above the line, and vice versa. All the models showing fragmentation 
are distributed above the solid line, while almost all the models showing jet are distributed 
below the solid line. The solid line clearly separates the models; fragmentation occurs but 
no jet appears when /3 > 7 , and jet appears after without fragmentation when (3 < 7 . In 
addition, the solid lines almost coincide with the boundary between the rotation-dominated 
model and the magnetically dominated models (see Table [1] for which category each model 
falls into). As a result, in the collapsing primordial cloud, the cloud evolution is mainly 
controlled by the centrifugal force than the Lorentz force when f3 > 70, while the Lorenz 
force is more dominant than the centrifugal force when 70 > (3q. 

The upper and right axes of Figure [10] mean the magnetic field strength and angular 
velocity of the initial cloud, respectively. They can be scalable at any initial density as 
(n c /10 3 cm -3 ) 2 / 3 , assuming the spherical collapse, which is approximated well in the early 
phase. A jet is driven when the initial cloud has magnetic field of B > 10 _9 (n c /10 3 cm~ 3 ) 2 / 3 G 
if the clou d rotates slowly as < 4 x 10~ 17 (?2 C /10 3 cm -3 ) 2 / 3 s" 1 . This condition corresponds 



to that in iMachida et al.l (j2006dl ) and is also consistent with the condition of equation ([I 



A strong jet is expected in primordial star formation. The power of a jet, e.g., a mass 
ejection rate, is considered to be controlled by the accretion rate as indicated in present-day 
star formation; the mass ejection rate of a jet is 1/10 of the mass accretion onto the protostar. 
The accretion rate of primordial star formation is expected to be considerably larger than 
that at present day, and it produces a stronger jet. The life time of the jet also seems to 
be controlled by accretion in the present-day; a jet stops when mass accr etion stops. For 



Population III stars, the gas accretion does not halt within their lifetimes (jOmukai fc Palla 



200ll . 120031 ) . Therefore, a jet also may continue during the all lifetime of the protostar, and 
the strong jet propagates to disturb a surrounding medium significantly. The disturbance of 
the medium could trigger the subsequent star formation as frequently observed in present-day 
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star formation. 



Assuming the power law growth of B oc n 2 J^ , the critical strength of the magnetic 
field, B = 1CT 9 G at n c = 10 3 cm~ 3 corresponds to B = 5 x 1CT 13 G at n c = 0.01 cm" 3 , 
which is much stronger than the background magnetic field of 10~ 18 G derived by llchiki et al. 
( 20061 ). However, when the magnetic field is amplified to B ~ 10 _9 (rz c /10 3 cm -3 ) 2 / 3 G by 



some mechanisms, the magnetic field can affect the collapse of the primordial cloud. Even if 
a cloud has a magnetic field weaker than the critical strengt h B n = 10~ 9 G, the magn etic field 



may play an important role after the protostar formation. iTan fc Blackmanl (120041 ) studies 
analytically the evolution of accretion disks around the first stars, suggesting that magnetic 
fields amplified in the circumstellar disk eventually give rise to protostellar jets during the 
protostellar accretion phase. 

Rotation promotes fragmentation when the first collapsed objects has the angular veloc- 
ity of fio ^ 10~ 17 (n c /10 3 cm -3 ) 2 / 3 s _1 as shown in Figure [TOl This condition coincides with 
that given by equation (|T6l) . The fragmentation is expected to produces binary or multiple 
stellar system. When a multiple stellar system is formed, some stars can be ejected by close 
encounters. At the protostar formation epoch, the protostar has a mass of M ~ 1O~ 3 M . 
The ejected proto-Population III stars may evolve to metal-free brown dwarfs or low-mass 
stars. When a binary component in a multiple stellar system is ejected from the parent 
cloud by proto stellar interaction, a low-mass metal free binary may also be appeared in the 
early universe. ISuda et al.l (120041) indicat es that the extremely metal-poor ([Fe/H]<-5) stars 



(IChristlieb et al 



2001 



Frebel et al.ll2005l ) discovered recently are formed as binary members 



from metal-fr ee gas, and then hav e been polluted by the companion stars during the stel- 
lar evolution. iKomiya et al.l (120061 ) shows that a binary frequency in Population III star is 
comparable to or larger than that at present day. In order to confirm the ejection scenario, 
the further long-term calculations are required. 
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Fig. 1. — (a-/) Time sequence of model 23 [(/3 , 70) = (1CT 3 , 2 x 1CT 3 )] before the pro- 
tostar formation n c < 10 21 cm" 3 . In each panel, the structures of the high-density region 
(red isosurface) and magnetic field lines (black and white streamlines) are plotted in three 
dimensions, while the density contours (false color and contour lines) and velocity vectors 
(thin arrows) are projected in each wall surface. The central number density n c , elapsed 
time t, grid level I, and grid size are also shown in each panel. 
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Fig. 2.— Time sequence of model 23 [(0 O , 70) = (1(T 3 , 2 x 10 -3 ) ] with I = 28 level of grid 
after the protostar formation n c > 10 21 cm -3 . The structure around the center of the cloud 
is plotted on z = plane (upper panels), and y = plane (middle panels), and in three 
dimension (lower panels). The density distribution (false color and contours) and velocity 
vectors (arrows) are plotted in the upper and middle panels. The red curves denote a border 
between inflow and outflow regions in the upper panels. The magnetic field lines (black and 
white streamlines) , high-density region (red isosurface), and outflow region (transparent red 
isosurface) are plotted in each lower panel. The density contours (false color and contour 
lines) and velocity vectors (thin arrows) are also projected on each wall surface in each lower 
panel. The central number density n c , and elapsed time t are denoted at the top of each 
upper panel. The elapsed time after the protostar formation t c is shown at the bottom of 
each upper panel. 
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Fig. 3. — Time sequence of model 21 [(/3 , 7o) = (10~ 4 , 2 x 10~ 7 ) ]. The density distribution 
(false color and contours) and velocity vectors (arrows) on the z = plane are plotted in 
each upper panel. The magnetic field lines (black and white streamlines) and high-density 
region (red isosurface) in three dimension are plotted in each lower panel. In each lower 
panel, the density contours (false color and contour lines) and velocity vectors (thin arrows) 
are also projected on each wall surface. The central number density n c , elapsed time t, and 
elapsed time after the protostar formation t c are denoted in each upper panel. The grid level 
is shown in the upper left corner of each panel. 
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Fig. 4. — Final states on the z = plane against parameters 70 and (3q. The model numbers 
are described in the upper left corner outside each panel. The density distribution (color- 
scale) is plotted in each panel. The grid level, maximum number density (n), and grid scale 
are denoted inside each panel. 70 - /3 plane is divided into four regions indicated by colors: 
the fragmentation (blue), non-fragmentation (red), merger (green), and no-collapse (gray) 
models. 
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Fig. 5.— Final states on the y = plane against parameters 70 and (3q. The model 
numbers are described in the upper left corner outside each panel. Density {false color and 
contours), velocity vectors (arrows) are plotted in each panel. The red line represents the 
border between the infall and outflow. The central number density n c , elapsed time t, grid 
level, grid scale and velocity unit are shown in each panel. 
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Fig. 6. — Final states in three dimensions against parameters 70 and (3q. The model numbers 
are described in the upper left corner outside each panel. The magnetic field lines (black and 
white streamlines) , high-density regions (red isosurface) , and outflow region (transparent red 
isosurface) are plotted in each panel. The density contours (false color and contour lines) 
and velocity vectors (thin arrows) are also projected in each wall surface. The grid level, 
and grid scale are shown in each panel. 
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Fig. 7. — Magnetic flux-Spin relation. The evolution of magnetic flux density and angular 
velocity at the cloud center are plotted. The x-axis indicates the square root of the magnetic 
pressure [B c /{M) l/2 ] divided by the square root of the thermal pressure (P c ). The y-axis 
represents the angular speed (f2 c ) divided by the free-fall rate [(AirGpc) 1 ^ 2 ]. The symbols 
represent the magnetic field and the angular velocity at the initial state (*) and at each central 
density (o: n c = 10 5 , 10 8 , 10 11 • • cm -3 ). Each line denotes the evolution path from the initial 
state (n c = 1.86 x 10 3 cm -3 ). The arrows mean the direction of evolution. The thick grey 
band denotes the magnetic flux-spin relation f^/[(0.25) 2 x 47rGp c ] + B 2 /[(0.2) 2 x 8ttP c ] = 1 
[see eq. ffTD ]. 
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Fig. 8. — Angular velocity (Q c ) normalized by the freefall timescale (4irGp c ) against the 
central number density n c for models 3, 9, 15, 21, 27, and 33. The relation £l c / '(AirGpc) 1 ^ 2 oc 
nH % is also plotted. 
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Fig. 9.— The magnetic field strength normalized by the square root of the density 
[B z /(8irc 2 gPc) 1 / 2 ] (upper panel), where c S) o is the initial sound speed, and inverse root of 
the plasma beta /3p~ 1//2 (lower panel) are plotted against the central number density n c for 
models 31-35. The relations B c /(4irc 2 p c ) 1 ^ 2 oc n l J % (upper panel) and f3p 1 ^ 2 oc n° c 12 (lower 
panel) are also shown. 
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Fig. 10. — Fragmentation and jet conditions for parameters 70 and (3q. Symbols mean 
the models showing fragmentation (o), jet (□), both fragmentation and jet (A), neither 
fragmentation nor jet (x), and no collapse (o). The upper and right axes indicate the 
initial magnetic field strength B , and angular velocity f2 depending on the initial density 
(n Cj o/10 3 cm _3 ) 2//3 , respectively. The thick line represent (3q = 70. Bottom axis means the 
mass to magnetic flux ratio (M/$) norm normalized by its critical value. 



